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We investigate the dynamics of closed FRW universe and anisotropic Bianchi type-IX universe 
characterized by two scale factors in a gravity theory including a higher curvature (Gauss-Bonnet) 
term. The presence of the cosmological constant creates a critical point of saddle type in the phase 
space of the system. An orbit starting from a neighborhood of the separatrix will evolve toward the 
critical point, and it eventually either expands to the de Sitter space or collapses to the big crunch. 
In the closed FRW model, the dynamics is reduced to hyperbolic motions in the two-dimensional 
center manifold, and the system is not chaotic. In the anisotropic model, anisotropy introduces 
the rotational mode, which interacts with the hyperbolic mode to present a cylindrical structure of 
unstable periodic orbits in the neighborhood of the critical point. Due to the non-integrability of 
the system, the interaction of rotational and hyperbolic modes makes the system chaotic, making 
it impossible for us to predict the final fate of the universe. We find that the chaotic dynamics 
arises from the fact that orbits with even small perturbations around the separatrix oscillate in the 
neighborhood of the critical point before finally expanding or collapsing. The chaotic character is 
also evidenced by the fractal structures in the basins of attraction. 

PACS numbers: 98.80.Jk, 04.50.Kd, 05.45.Pq 



I. INTRODUCTION 

The pioneering work of Belinskii, Khalatnikov, and Lif- 
shitz [T] has shown that anisotropic Bianchi type-IX uni- 
verses could exhibit chaotic dynamics as it approaches 
the cosmological singularity (£ — » 0). The anisotropic 
Bianchi IX universe, as it approaches the singularity, ex- 
hibits an oscillatory mode consisting of an infinite se- 
quence of Kasner eras, during which two scale factors os- 
cillate and the remaining scale factor decreases monoton- 
ically. Independently, Misner [2] also suggested a chaotic 
approach to initial singularity in his Mixmaster universe 
(vacuum Bianchi IX universe characterized by three scale 
factors) . 

While the aspect of chaotic dynamics helped deepen 
our understanding of singularities in general relativity, 
there have been significant challenges in finding an invari- 
ant measure of chaos. A standard measure of chaos is the 
Lyapunov exponent, and a debate started on whether or 
not the Mixmaster universe is chaotic, when some studies 
found zero Lyapunov exponent for the Mixmaster uni- 
verse, while others found a positive Lyapunov exponent 
[3H5]. The conflict was resolved when it was realized 
[B] that Lyapunov exponent is coordinate dependent and 
therefore is not a reliable measure of chaos in general 
relativity. 

Cornish and Levin [7] demonstrated the appearance 
of chaos in the dynamics of the Friedmann-Robertson- 
Walker (FRW) models with a cosmological constant and 
scalar fields conformally coupled to the geometry. The 
presence of the cosmological constant creates a saddle 
point in the phase space, and the separatrix connects 
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this saddle point to other critical points. When there are 
interactions between the geometry with the scalar fields, 
the separatrix breaks up and becomes fractal. The work 
by de Oliveira, Soares, and Stuchi [5] emphasized that 
even small perturbations induces the breaking of highly 
unstable separatrix. The presence of a positive cosmolog- 
ical constant and a perfect fluid creates a saddle-center 
in the phase space, and non-integrability of the system 
induces distortion and twisting of the topology of homo- 
clinic cylinders in the neighborhood of the critical point. 
The fractal and topological structures are coordinate- 
independent and thus provide invariant characterization 
of chaos in relativistic theories. 

Regarding such a system as a model of the early uni- 
verse, a natural question that arises is whether the Ein- 
stein gravity is reliable or not. It has been pointed out 
[3] that we might have the trans-Planckian issue; at early 
stages of inflation, the quantum gravity effects are not 
necessarily negligible. Previous works on chaotic dynam- 
ics in the per-inflationary era focused on Einstein gravity, 
but as we approach the Planck scale, it is reasonable to 
expect higher order corrections to the Einstein-Hilbert 
action. Superstring theory is the leading candidate for 
a description of physics at such scale. Our goal in this 
paper is to use an invariant characterization of chaos and 
investigate the existence of chaotic dynamics in "string- 
inspired" modified gravity. Our model takes into account 
the effects of the higher curvature terms that typically 
arise in the one-loop low-energy effective superstring ac- 
tion. For simplicity we focus on a model in which a scalar 
field representing the string moduli is non-minimally cou- 
pled to the Gauss-Bonnet curvature. For a related work 
on chaotic dynamics of higher curvature modified gravity, 
see e.g. [TO] , 

The Bianchi type-IX models with a scalar field coupled 
to the Gauss-Bonnet term are studied in somewhat dif- 
ferent context of singularity avoidance in . The new 



features introduced by our model are a positive cosmo- 
logical constant and a perfect fluid, which create saddle 
points in the phase space of the system. Because this 
method does not work for the case of zero cosmological 
constant, our work implies that in order for the universe 
to inflate we need to include a positive cosmological con- 
stant in the low-energy effective action. We introduce 
this modified action in Sec. [TTJ We first consider the 
closed FRW universe in Sec. where the degrees of 
freedom are the scale factor and a scalar field, and dis- 
cuss the basic characteristics of the closed FRW model. 
We introduce anisotropy in the metric in Sec. |IV| and 
discuss how anisotropy creates a topological structure of 
cylinders near the critical point. In Sec. |V| we present 
numerical evidence of cylindrical topology, oscillatory be- 
havior around the critical point, and fractal structures 
in the basins of attraction when small perturbations are 
introduced in the metric and/or the scalar field. The 
fractal structures in the basins of attraction, as well as 
the topology of cylinders, constitute invariant character- 
ization of chaos, and we conclude that chaotic dynamics 
can arise in our string-inspired model. 



II. ONE-LOOP EFFECTIVE ACTION 



We start with the action in the Einstein frame given 



by HMS] 



S = 



d A x, 



^matter (1) 



where R and a are the Ricci scalar curvature and a scalar 
field, respectively. In our units, the gravitational con- 
stant corresponds to G — l/8ir. A positive cosmological 
constant is necessary for the existence of a saddle point 
in phase space and the universe to inflate. The Gauss- 
Bonnet curvature is given by 

R-gb — R 2 — 4R al3 R a i3 + R a ^ 1& 'R a p 7 s (2) 

and the function £ (cr) depends on details of string theory 
compactified geometry jTT]. For example, when cr is the 
dilaton, we have 



£(cr)=e CT 



(3) 



In type II superstring, a can be regarded as the mod- 
ulus of compactified dimensions. The form of £ (cr) in 
a particular compactification can be found in |12l ITT] . 
The function £ (cr) determines the coupling of cr and the 
geometry, and is expressed with Dedekind rj function as 



£ (a) = - In [2eV (i 
Due to the modular property 



(4) 



(5) 



the function £ (cr) is even in cr, has a global minimum 
at a = 0, and increases exponentially as a — ► ±00. We 
will be interested in the small a behavior in the follow- 
ing sections. We assume for simplicity that £ (cr) can be 
approximated as 



UP) 



1 



(6) 



III. CLOSED 
FRIEDMANN-ROBERTSON- WALKER 
UNIVERSE 

We first consider a cosmological model characterized 
by the scale factor a (i) with the line element given by 



ds 2 



dt 2 + a 2 (t) {u 1 ) 2 + {u> 2 ) 2 + (u 3 Y 



(7) 



The Bianchi type-IX invariant one-forms uj 1 are given by 
uj 1 = sin ip d9 — sin 9 cos ip d(f> 



oj 2 = cos ifi d9 + sin 9 sin ip d(f> 
oj 3 = cos 9d(f> + dip 



(8) 



which are chosen to satisfy du; 1 = — j^uji A ui k and 
e l jk is the completely antisymmetric tensor (e 1 23 = l)- 
This represents the closed FRW universe with positive 
curvature k = +1. We assume that the matter content 
is a perfect fluid. The energy- momentum tensor of the 
perfect fluid can be written in the form 



(P + P) $0^0 + PfJ 



(9) 



where p and p are energy density and pressure, respec- 
tively. For simplicity, we assume that the perfect fluid 
may be represented by "dust," that is, 7 = in the equa- 
tion of state p = j p. However, even for general perfect 
fluids, we can expect similar features [181 . 
For the action ([lj the Lagrangian is given by 



L = -a — Aa 3 — 3aa 2 
4 



1 



a cr 



A 



+ -crc>(3d + 4a 3 ) (10) 
8 



where the overdot denotes differentiation with respect to 
time t. The equations of motion are 



Aa 2 
1 



1 



a 2 a 2 — 2aa 



A (a 2 



4d 2 cr 2 



8crcrdd + aa + Aaah z ) = 



3a 2 dcr + a 3 cr + -Acrd (1 + 4d 2 ) =0 (11) 
8 

with the Hamiltonian constraint given by 

3 3 1 

H = -3aa 2 --a+a 3 A+-\aaa (l + 4:a 2 ) + -a 3 a 2 +E = 

4 8 2 

(12) 



where Eq = pa 3 corresponds to the total matter content 
of the model. 

We can always decompose a set of second-order differ- 
ential equations into a set of first-order differential equa- 
tions by redefining variables. Thus, we get a set of four 
coupled first-order ordinary differential equations 



dt 



a = a, 



dt 



a = a 



dt 



a = 8a 2 



3Actcto, (1 + 4d 2 



dt 



+ a 3 (8A - Air 2 ) + a (l + 4d 2 ) (-2 + Act 2 ) 



<? = \ - 384a 3 ctCT + 3Act ( - (l + 4a 2 ) 2 (-2 + Act 



4a 2 [-2A + ct 2 + 4a 2 (-2A + 5ct 2 )] ) U (13) 



where x is given by 

X 



1 



64a 3 (2a - Xaaa) + 3A 2 ct 2 (1 + 4a 2 Y 

) admits the stati 
i is just the critical point P. 



The dynamical system ( 13 1 admits the static Einstein 
universe as a solution, wl 



Its coordinates are given by 



P : a — ao = 



1 



, ct = ctq , a = 0, ct = 



2VA 

where cto is a constant. The critical energy is given by 

1 



(15) 



E r . 



4vA 



(16) 



When the scalar field is absent and does not interact 
with the curvature, the equations of motion can be in- 
tegrated exactly, and the dynamical system admits an 
invariant manifold Ai. The invariant manifold M is de- 
fined by 



M : ct = 0, ct = 



(17) 



and the dynamical system simplifies to an exactly inte- 
grable two-dimensional system 



dt 



a = a 



dt a 8a 2a 



-Aa 



(18) 



In Fig. [T]we display the phase space portrait of the in- 
variant manifold A4 in the (a, a) plane. We introduce the 
conformal time dr/ = dt/a so that the dynamical system 
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Figure 1: Phase space portrait of the invariant manifold M 
in conformal time dr; = At/a. The orbit S is the separatrix 
and the point P is the critical point. 



(14) (18 1 becomes regular at a = 0. The phase space is di- 



vided by a separatrix into two types of orbits: those that 
collapse into the big crunch and those that expand into 
the de Sitter space. The critical point P intersects the 
invariant manifold M at cto — 0. Furthermore, there are 
attractors corresponding to stable and unstable de Sitter 
spaces. When the dynamics is dominated by the cosmo- 
logical constant A, it can be shown that the scale factor 
a (t) expands exponentially and approaches the stable de 

Sitter attractor as a (t) ~ eV /A / 3 * . 



Let us now linearize the dynamical equations (13 1 
about the critical point P. We move the critical point 
to the origin by redefining 



a + ao, ct— >ct + 
+ ct^ct + 



then we obtain 




/higher^ 
I order 
terms 



(19) 



(20) 



where the constant matrix associated with linearizing the 



system ( 13 ) about the critical point P is given by 



An 



/0 1 0\ 

1 

A 

\0 0/ 



The matrix (211 has four eigenvalues 

Ai 2 = ±VI, A 3 , 4 = 2 



(21) 



(22) 



According to the center manifold theorem [TJ, the 
study of the dynamics near the critical point can be re- 
duced to the study of the dynamics restricted to the asso- 
ciated two-dimensional invariant manifolds W c (0) near 
the critical point. Without actually calculating the cen- 
ter manifold, we argue as follows. Under the coordinate 
transformation 



(23) 



(° 





1 






1 











u 2 
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1 


1 "3 


Vo 


1 





o J 


\u 4 



the matrix (21) assumes the Jordan canonical form 



J 



o 








-Va 











(24) 



Thus the four-dimensional system ( 13 I restricted on the 
center manifold W c (0) is a two-dimensional vector field 
whose linear part is given by 



1 





(25) 



Bogdanov [20] has shown that the system which has a 
linear part (25) is locally topologically equivalent near 



the critical point to the normal form 



v 2 



V2 = 01 + P2V2 + Vx + JV1V2 



(26) 



where 7 = ±1. The dynamics of (26 1 is not qualita- 
tively changed by the higher order terms in the normal 
form. The normal form ( 26 1 has no homoclinic or peri- 



odic solutions for j3\ > 0. If we rescale the variables and 
parameters as 



vi = e 2 qi, v 2 = e 6 pi, Pi 

where e > and rescale the time as 

t 

t - 



4 

e /I, 



(27) 



(28) 



where fi is a parameter with /1 > corresponding to the 
near-zero eigenvalues being real. For /i < there exist 
homoclinic solutions. However, we restrict our attention 
to only the hyperbolic motions and do not consider the 
homoclinic orbits, because homoclinic orbits extend to 
a < region of the phase space, which does not have a 
physical meaning. 

In the parlance of nonlinear dynamical systems theory, 
this is the Bogdanov- Takens bifurcation. To perform a 
global analysis that includes the effect of the O (e) part 
of (29 1 on this integrable structure, Melinkov's method 
[H?] can be used. A detailed description of bifurcation 
diagrams and phase portraits of the Bogdanov- Takens 
bifurcation can be found elsewhere (cf. [21 ) and need 
not be repeated here. 

What the above linear analysis shows is that although 
we could expect simple hyperbolic motion of the planar 
Bogdanov- Takens bifurcation, we do not expect chaotic 
dynamics in our closed FRW model. As we will see in 
Sec. 



IV 



the chaotic dynamics, especially the rotational 
mode of the oscillatory Mixmaster dynamics, arises from 
the presence of imaginary eigenvalues of the linearized 
matrix. 



IV. AXISYMMETRIC BIANCHI IX UNIVERSE 

Chaos is expected to come from the anisotropy of the 
universe. For simplicity, we restrict our attention to two 
distinct scale factors a (t) and b (t) and consider an ax- 
isymmetric Bianchi IX cosmological model with the line 
element 

ds 2 = -dt 2 + a 2 (t) [(uj 1 ) 2 + (co 2 ) 2 ] + b 2 (t) (uj 3 ) 2 (31) 
For the action the Lagrangian is given by 



L = b- Aa 2 b - - ba 2 - 2aab + ^-a 2 b& 2 



4a 2 
A . / b 3 d 



1 

2' 
.tfb 
1~ 2 



ad I 2— +46- 3^ +Aa 2 b\ (32) 



and the equations of motion are 



then the system becomes 



<7i =Pi 



Pi = M + Qi + e (Pi + 7<?iPi) 



(29) 



b 3 

2Aab - 2db - aba 2 - 2ab - 2ab 

2a fi 



+ A 



4a3 



+ db J (a 2 + cer) + aa (db + db 



= 



For e = 0, the rescaled equations become exactly inte- 
grable Hamiltonian system with the Hamiltonian given 

by 



H (qi,pi) 



Pi 



Q 3 i 



Mi 



(30) 



l+Aa 
+ A 



36^ 
" 4a 2 
1 d 2 
2 + Y 



2 

3^ 

8a 2 



2 • 2 
a a 



2ad 



= 



2abaa + a 2 b& + a 2 ba 



the matrix ( 37 1 assumes the Jordan canonical form 



b 3 a 2 b 2 ab bb 2 b 3 & 
i— r- + 12-t- - 6-^- + 2^- 



8a6a + 46-3- 



b 2 b 



Aa 2 b 



(33) 



with the Hamiltonian constraint given by 

4a 2 
b 2 b 



H 



-a b — 2adb — b 



a 2 bA 



1 



a 2 ba 2 



X ( ab 3 



46-3- 



12a^6 + £ = (34) 



where _E = pa 2 6 corresponds to the total matter content 
of the model. 

We can always decompose a set of second-order differ- 
ential equations into a set of first-order differential equa- 
tions by redefining variables. Thus we get a set of six 
coupled first-order differential equations in the variables 

I a, b, cr,a,b,&). The dynamical system (33 1 admits the 



static Einstein universe as a solution, which is just the 
critical point P. Its coordinates are given by 

P : a = b = ao = / _ , a = <tq, a = b = & = (35) 



2VA 

where er is a constant. The critical energy is given by 

1 



4VA 



(36) 



When the scalar field is absent and the isotropy is re- 
stored, the equations of motion are exactly integrable and 
the dynamical system admits the same invariant manifold 
A4 of Fig. [I] As before, the critical point P intersects 
the invariant manifold Ai at cjq = 0. 

The constant matrix associated with linearizing the 



system ( 33 1 about the critical point P is given by 
/ 



An 



Q 



-2A 
6A 




1 0\ 

1 

1 

3A 

-5A 

0/ 



(37) 



The matrix ( 37 1 has six eigenvalues 
Ai o = ±VA, A 3 , 4 = 2 , 



A, 



±2iV2A 



(38) 



Compared to the four eigenvalues (22 1 in the isotropic 
case, the anisotropy in the metric has produced an ad- 
ditional pair of imaginary eigenvalues. Without loss of 
generality, we fix A = 1/4 so that ao = 1. Under the 
transformation 
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o / 



(40) 



Thus, with canonical transformation, the Hamiltonian to 
the quadratic order can be expressed as 



H ■ 



9 
Pi 



1 



(pi 



4 



Eq — E cr it 
+ 0(3) = O (41) 



In a small neighborhood of the critical point, the higher 
order terms in the expansion O (3) is negligible and we 
may assume that the energy Eq — E cr a is small. Then 
the Hamiltonian may be approximated as 



H 



1 



'2 2\ i 2 
J>2 - 12) +Pl 



1 

71 



(pi 



?3 2 ) 



In this linear regime, the Hamiltonian (42) is separable. 
If we define the partial energies as 



Eq — E cr it — 

(42) 



Ei=p{ 



E hyp = ^ (pt - ql) 



Ey 



1 

71 



>3 + ft") 



they are approximately conserved separately, 

E x 



E h 



Eq 



(43) 



(44) 



(45) 



(46) 



The topological structure near the critical point cre- 
ated by these separable partial energies in the case of gen- 
eral relativity was first described by de Oliveira, Soares, 
and Stuchi in [S\. Compared to their model, we have 
an additional partial energy E±. Ignoring the energy 
Ei for the moment, let us concentrate on the hyper- 
bolic motion energy Ei iyp and the rotational motion en- 
ergy E rot . If Ehyp — 0, we have either P2 — 02 = or 
P2 = ±92 • When p 2 = 02 = 0, we have p 2 + q 2 = constant 
and the motion will be described by periodic orbits te 
in the (pi,qi) plane. These periodic orbits depend on 
the parameter Eq. When p 2 = ±02, the motion will 
be described by linear stable V s and unstable V u one- 
dimensional manifolds in the (^2,92) plane. In addition, 
we have E\ = p\ = constant, which generates a linear 
one-dimensional manifold T e ■ The direct product of te 
with V s , V u , and r^ generates the topological structure 
of stable cylinders te X Te xV s and unstable cylinders 
t e x r_e x Kt- Thus the flow in the phase space will be 
E. 2 x S 1 . 



The center manifold is the nonlinear extension of the 
linear regime that corresponds to Eh yp = 0. The inter- 
section of this center manifold with the energy surface 
Eq = E cr it corresponds to p\ — p2 = 92 = 0, which is 
just the critical point P. Since E rot and E\ is always 
positive, for Eq > E cr n the center manifold does not in- 
tersect the energy surface in the linear regime. However, 
when we include nonlinear terms, the Gauss-Bonnet term 
may have positive as well as negative energy. Therefore, 
contrary to the general relativity case, the center mani- 
fold can intersect the energy surface even if E Q > E critl 
and we may have a saddle structure for Eq > E cr u . 

In general, an orbit that approaches the neighborhood 
of the critical point P will have E\ ^ 0, Eh yp 7^ 0, and 
E ro t 7^ 0. As this orbit with energy E cr i t — Eo approaches 
the critical point, the non-integrability of the system ( 33 ) 



makes it impossible to predict the amount of energy that 
will be partitioned into each mode E±, Eh yp , or E rot . 
Since we cannot predict how much energy will be trans- 
formed into which mode, we cannot predict whether the 
orbit will collapse into the big crunch or escape to infla- 
tion. The rotational mode arising from anisotropy corre- 
sponds to the oscillatory mode of the Mixmaster dynam- 
ics and is crucial in this non-predictability. 

In the parlance of nonlinear dynamical systems theory, 
the double zero eigenvalues lead to Bogdanov-Takens bi- 
furcation and a pair of purely complex eigenvalues lead to 
Andronov-Hopf bifurcation [TU] in the four-dimensional 
center manifold. We have new possibilities in our system 
due to the interaction of the Bogdanov-Takens bifurca- 
tion with the Andronov-Hopf bifurcation. The simple hy- 
perbolic motion of the planar Bogdanov-Takens bifurca- 
tion transversally intersects two-dimensional stable and 
unstable manifolds of periodic orbits, and this interaction 
leads to chaotic dynamics. 



V. NUMERICAL RESULTS 

In our numerical experiments, we closely follow the 
approach of de Oliveria, Soares, and Stuchi [5]. We fix 
A = 0.25 so that the coordinates of the critical point 
P are given by a = b = 1.0, er = 0, a = b = a = 0, 
and E cr it = 0.5. Let Sq be a point on the separatrix 
S in the invariant manifold Ai (a = & = 0). For the 
results presented in this paper, we choose the coordi- 
nates a = b = 0.9 and a = b = 0.051818 772 5, but 
any point near the separatrix would yield similar results. 
Around So, we perturb the separatrix in five variables 
(three in the closed FRW case) by an arbitrarily small 
amount S = 10~ 4 and use the Hamiltonian constraint to 
fix the remaining variable. The energy of the orbit Eq 
is chosen to be very close to the energy of the separa- 
trix so that the difference in energy AEq = E cr i t — Eq 
is much smaller than the perturbation 5. This mim- 
ics the uncertainty in the initial conditions. In physical 
terms, these initial conditions represent small perturba- 
tions in the scale factor and/or the scalar field in the 




Figure 2: Escape to inflation of 100 orbits around a point 
on the separatrix with coordinates a — b — 0.9, a — b = 
0.051818 772 5, and a = a = 0. The coupling constant is 
A = 16. The energy surface is given by E Q = 0.499 999 999 9 
and the radius of the sphere of initial conditions is S = 10~ 4 . 
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Figure 3: Collapse of 100 orbits around the same point in Fig. 
[|with A = 16, E Q = 0.499 999 980 2, and S = 10" 4 . 



pre-inflationary era. 

When the orbits are evolved starting from the point 
Sq, the orbits evolve toward the critical point since the 
initial conditions were chosen near the separatrix. After 
passing through the critical point, we expect the orbits 
to either collapse into the big crunch or expand to de 
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Figure 4: A magnified view of 100 orbits around the crit- 
ical point in the (a, a) plane for S = 10 -4 and Eo = 
0.499 999 997 5. 



Sitter space. The two possible outcomes for 100 orbits 
with S = 10~ 4 are shown for the inflation in Fig. [2] 
and for the collapse in Fig. [3j The final fate of the 
orbits depends on the energy E , and there exists an 
upper bound on Eq for which all orbits collapse and a 
lower bound on Eq for which all orbits escape to inflation. 
For an energy between this interval, some orbits collapse 
and other orbits inflate, resulting in an indeterminate 
outcome. 

In Fig. [4] we display a magnified view of the region 
around the critical point in the (a, d) plane. Note that 
the orbits oscillate around the separatrix as well as the 
critical point. This oscillatory mode is crucial in the mix- 
ing of the boundaries and the existence of chaotic dynam- 
ics. As this orbit with energy E cr n — Eq approaches the 
critical point, the non-integrability of the system makes 
it impossible to predict the amount of energy that will be 
partitioned into each mode Ei, Eh yp , or E rot . Thus, the 
outcome has a sensitive dependence on initial conditions. 

We now proceed to find a set of initial values that lead 
to an orbit approaching the de Sitter attractor. Such a 
set is called a basin of attraction. Here we use the stan- 
dard method and do a pixel-by-pixel computation of a 
400 x 400 grid. Since we are dealing with basins embed- 
ded in a six-dimensional phase space, we are forced to 
consider lower dimensional slices, and here we choose the 
(a, &) plane although similar fractal basins of attraction 
can be obtained for other slices in the phase space. As 
shown in Fig. [5] in the closed FRW universe, we have a 
sharply divided separatrix even if we have a contribution 
from the Gauss-Bonnet term. When the metric becomes 
anisotropic, the basins of attraction become highly frac- 
tal as shown in Fig. p\ A magnification of the inner 
region is shown in Fig. u\ and reveals self-similar fractal 





cr 



Figure 5: The basins of attraction in the (a, a) plane for the 
closed FRW model. The slice is through a = 0.9 and a is fixed 
by the Hamiltonian constraint. We chose the parameters A = 
16 and E = 0.499 999 997 5. The black regions correspond to 
collapse and the white regions inflation. 
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Figure 6: The basins of attraction in the (a, a) plane for the 
axisymmetric Bianchi IX model. The slice is through a — b = 
0.9, a = 0.051818 772 5, and b is fixed by the Hamiltonian 
constraint. The other parameters are the same as in Fig. [5] 



structure. Note that the black regions correspond to or- 
bits that collapse to the big crunch and the white regions 
to orbits that expand to the de Sitter space. 

The existence of fractal structures is not restricted to 
one particular plane and can be seen in other slices in 
the six-dimensional phase space. Fig. [8] shows that simi- 
lar fractal structures can also be seen in the (a, d) plane 
for universes with parameters and conditions identical to 
those used in Fig. [6] These fractal structures in the phase 
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Figure 7: A magnified view of Fig. [6] 
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Figure 8: The basins of attraction in the (a, a) plane. 



space cannot be removed by coordinate transformation 
and thus provide invariant characterization of chaos in 
our model. 



description of pre-inflationary stages of the Universe af- 
ter the Planck era. Higher curvature terms arise as the 
next-to leading terms in the superstring effective action, 
or from the renormalization of the stress tensor. We 
have included such effects in the form of a Gauss-Bonnet 
curvature term non-minimally coupled to a scalar field. 
The main features of our model are a positive cosmolog- 
ical constant and a perfect fluid, which produce a sad- 
dle point in the phase space. Due to the presence of 
the saddle point, an orbit starting in the neighborhood 
of the separatrix has two asymptotic possibilities, the 
de Sitter inflation and the big crunch collapse. In the 
closed FRW model, the dynamics near this critical point 
is reduced to simple two-dimensional hyperbolic motion 
of the Bogdanov-Takens bifurcation. Consequently, the 
closed FRW model is not chaotic. In the anisotropic 
Bianchi type-IX model, we restricted our attention to 
the axisymmetric case. The introduction of anisotropy 
produces unstable periodic orbits, which interact with 
simple hyperbolic motion of the planar Bogdanov-Takens 
bifurcation to produce chaotic dynamics. 

Our results extend the work of de Oliveira, Soares, and 
Stuchi in [5], who describe chaotic exit to inflation for the 
axisymmetric Bianchi IX universe in general relativity. 
In subsequent papers fTSJ [221 EI] the authors discuss the 
existence of "homoclinic chaos" in general relativity. Our 
model does not have such homoclinic orbits because of 
nonlinearity due to the coupling of Gauss-Bonnet curva- 
ture and the scalar field. Furthermore, extending a cos- 
mological model to the a < region of the phase space 
is physically unreasonable as pointed out in |24| . 

In the model considered in [5] , it is possible to find an 
energy E$ such that a small perturbation to a point on 
the separatrix is a chaotic set. As we have shown in this 
paper, this is also true when we include a stringy cor- 
rection to the Einstein-Hilbert action. We have shown 
that a small fluctuation in initial conditions leads to in- 
determinate outcome between collapse or inflation. Fur- 
thermore, we found some numerical evidence of fractal 
structures in the basins of attraction. The fractal basins 
of attraction, together with topology of cylinders near 
the critical point, are invariant characterization of chaos 
in our model. 
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